[Position Effect Variegation: Role of the Local Chromatin Context in Gene Expression Regulation].

Position effect variegation (PEV) is a phenomenon wherein the expression level of a gene strongly depends on its genomic position. PEV can be observed when a gene is moved via a chromosome rearrangement or identical genetic constructs are inserted into different regions of the genome. The eukaryotic genome has a domain organization, and gene activity within a domain depends not only on the nucleotide sequence of a gene, but also on the state of surrounding chromatin, thus being regulated epigenetically. Chromatin is a complex of DNA, RNA, and associated structural and regulatory proteins. The epigenetic status of chromatin depends on the replication time of a given genomic region, particular regulatory DNA motifs, and contacts with the inner nuclear envelope (lamina) and other chromosome regions (topologically associated domains). PEV results from the changes in the epigenetic state of a gene and provides a unique tool to study the molecular and biochemical processes that underlie the establishment and switching of epigenetic states. Understanding the molecular mechanisms of PEV in human is of clinical importance, in particular, for the detection and treatment of retroviral infections because the local chromatin state may determine the latent/active state transition of an infection, such as HIV. In addition, a large number of human neurodegenerative diseases are caused by epigenetic gene inactivation due to expansion of short repeats. Finally, to apply gene therapy methods, it is important to develop approaches that ensure a necessary level of transgene expression with sufficient accuracy.

Transcriptome Analysis in Yeast Reveals the Externality of Position Effects.

The activity of a gene newly integrated into a chromosome depends on the genomic context of the integration site. This "position effect" has been widely reported, although the other side of the coin, that is, how integration affects the local chromosomal environment, has remained largely unexplored, as have the mechanism and phenotypic consequences of this "externality" of the position effect. Here, we examined the transcriptome profiles of approximately 250 Saccharomyces cerevisiae strains, each with GFP integrated into a different locus of the wild-type strain. We found that in genomic regions enriched in essential genes, GFP expression tended to be lower, and the genes near the integration site tended to show greater expression reduction. Further joint analysis with public genome-wide histone modification profiles indicated that this effect was associated with H3K4me2. More importantly, we found that changes in the expression of neighboring genes, but not GFP expression, significantly altered the cellular growth rate. As a result, genomic loci that showed high GFP expression immediately after integration were associated with growth disadvantages caused by elevated expression of neighboring genes, ultimately leading to a low total yield of GFP in the long run. Our results were consistent with competition for transcriptional resources among neighboring genes and revealed a previously unappreciated facet of position effects. This study highlights the impact of position effects on the fate of exogenous gene integration and has significant implications for biological engineering and the pathology of viral integration into the host genome.

Molecular effects of dADD1 misexpression in chromatin organization and transcription.

BACKGROUND: dADD1 and dXNP proteins are the orthologs in Drosophila melanogaster of the ADD and SNF2 domains, respectively, of the ATRX vertebrate's chromatin remodeler, they suppress position effect variegation phenotypes and participate in heterochromatin maintenance. RESULTS: We performed a search in human cancer databases and found that ATRX protein levels were elevated in more than 4.4% of the samples analyzed. Using the Drosophila model, we addressed the effects of over and under-expression of dADD1 proteins in polytene cells. Elevated levels of dADD1 in fly tissues caused different phenotypes, such as chromocenter disruption and loss of banding pattern at the chromosome arms. Analyses of the heterochromatin maintenance protein HP1a, the dXNP ATPase and the histone post-translational modification H3K9me3 revealed changes in their chromatin localization accompanied by mild transcriptional defects of genes embedded in heterochromatic regions. Furthermore, the expression of heterochromatin embedded genes in null dadd1 organisms is lower than in the wild-type conditions. CONCLUSION: These data indicate that dADD1 overexpression induces chromatin changes, probably affecting the stoichiometry of HP1a containing complexes that lead to transcriptional and architectural changes. Our results place dADD1 proteins as important players in the maintenance of chromatin architecture and heterochromatic gene expression.

Maximizing antibody production in a targeted integration host by optimization of subunit gene dosage and position.

Historically, therapeutic protein production in Chinese hamster ovary (CHO) cells has been accomplished by random integration (RI) of expression plasmids into the host cell genome. More recently, the development of targeted integration (TI) host cells has allowed for recombination of plasmid DNA into a predetermined genomic locus, eliminating one contributor to clone-to-clone variability. In this study, a TI host capable of simultaneously integrating two plasmids at the same genomic site was used to assess the effect of antibody heavy chain and light chain gene dosage on antibody productivity. Our results showed that increasing antibody gene copy number can increase specific productivity, but with diminishing returns as more antibody genes are added to the same TI locus. Random integration of additional antibody DNA copies in to a targeted integration cell line showed a further increase in specific productivity, suggesting that targeting additional genomic sites for gene integration may be beneficial. Additionally, the position of antibody genes in the two plasmids was observed to have a strong effect on antibody expression level. These findings shed light on vector design to maximize production of conventional antibodies or tune expression for proper assembly of complex or bispecific antibodies in a TI system.

The impact of genetic background and cell lineage on the level and pattern of gene expression in position effect variegation.

BACKGROUND: Chromatin-based transcriptional silencing is often described as a stochastic process, largely because of the mosaic expression observed in position effect variegation (PEV), where a euchromatic reporter gene is silenced in some cells as a consequence of juxtaposition with heterochromatin. High levels of variation in PEV phenotypes are commonly observed in reporter stocks. To ascertain whether background mutations are the major contributors to this variation, we asked how much of the variation is determined by genetic variants segregating in the population, examining both the level and pattern of expression using the fruit fly, Drosophila melanogaster, as the model. RESULTS: Using selective breeding of a fourth chromosome PEV reporter line, 39C-12, we isolated two inbred lines exhibiting contrasting degrees of variegation (A1: low expression, D1: high expression). Within each inbred population, remarkable similarity is observed in the degree of variegation: 90% of the variation between the two inbred lines in the degree of silencing can be explained by genotype. Further analyses suggest that this result reflects the combined effect of multiple independent trans-acting loci. While the initial observations are based on a PEV phenotype scored in the fly eye (hsp70-white reporter), similar degrees of silencing were observed using a beta-gal reporter scored across the whole fly. Further, the pattern of variegation becomes almost identical within each inbred line; significant pigment enrichment in the same quadrant of the eye was found for both A1 and D1 lines despite different degrees of expression. CONCLUSIONS: The results indicate that background genetic variants play the major role in determining the variable degrees of PEV commonly observed in laboratory stocks. Interestingly, not only does the degree of variegation become consistent in inbred lines, the patterns of variegation also appear similar. Combining these observations with the spreading model for local heterochromatin formation, we propose an augmented stochastic model to describe PEV in which the genetic background drives the overall level of silencing, working with the cell lineage-specific regulatory environment to determine the on/off probability at the reporter locus in each cell. This model acknowledges cell type-specific events in the context of broader genetic impacts on heterochromatin formation.

Monitoring of switches in heterochromatin-induced silencing shows incomplete establishment and developmental instabilities.

Position effect variegation (PEV) in Drosophila results from new juxtapositions of euchromatic and heterochromatic chromosomal regions, and manifests as striking bimodal patterns of gene expression. The semirandom patterns of PEV, reflecting clonal relationships between cells, have been interpreted as gene-expression states that are set in development and thereafter maintained without change through subsequent cell divisions. The rate of instability of PEV is almost entirely unexplored beyond the final expression of the modified gene; thus the origin of the expressivity and patterns of PEV remain unexplained. Many properties of PEV are not predicted from currently accepted biochemical and theoretical models. In this work we investigate the time at which expressivity of silencing is set, and find that it is determined before heterochromatin exists. We employ a mathematical simulation and a corroborating experimental approach to monitor switching (i.e., gains and losses of silencing) through development. In contrast to current views, we find that gene silencing is incompletely set early in embryogenesis, but nevertheless is repeatedly lost and gained in individual cells throughout development. Our data support an alternative to locus-specific "epigenetic" silencing at variegating gene promoters that more fully accounts for the final patterns of PEV.

Position effects on promoter activity in Escherichia coli and their consequences for antibiotic-resistance determinants.

The activity of any bacterial promoter is generally supposed to be set by its base sequence and the different transcription factors that bind in the local vicinity. Here, we review recent data indicating that the activity of the Escherichia coli lac operon promoter also depends upon its chromosomal location. Factors that affect promoter activity include the binding of nucleoid-associated proteins to neighbouring sequences, supercoiling and the activity of neighbouring promoters. We suggest that many bacterial promoters might be susceptible to similar position-dependent effects and we review recent data showing that the expression of mobile genes encoding antibiotic-resistance determinants is also location-dependent, both when carried on a bacterial chromosome or a conjugative plasmid.

Regulation of position effect variegation at pericentric heterochromatin by Drosophila Keap1-Nrf2 xenobiotic response factors.

The Keap1-Nrf2 signaling pathway plays a central role in the regulation of transcriptional responses to oxidative species and xenobiotic stimuli. The complete range of molecular mechanisms and biological functions of Keap1 and Nrf2 remain to be fully elucidated. To determine the potential roles of Keap1 and Nrf2 in chromatin architecture, we examined the effects of their Drosophila homologs (dKeap1 and CncC) on position effect variegation (PEV), which is a transcriptional reporter for heterochromatin formation and spreading. Loss of function mutations in cncC, dKeap1, and cncC/dKeap1 double mutants all suppressed the variegation of wm4 and SbV PEV alleles, indicating that reduction of CncC or dKeap1 causes a decrease of heterochromatic silencing at pericentric region. Depletion of CncC or dKeap1 in embryos reduced the level of the H3K9me2 heterochromatin marker, but had no effect on the transcription of the genes encoding Su(var)3-9 and HP1. These results support a potential role of dKeap1 and CncC in the establishment and/or maintenance of pericentric heterochromatin. Our study provides preliminary evidence for a novel epigenetic function of Keap1-Nrf2 oxidative/xenobiotic response factors in chromatin remodeling.

Using a phiC31 "Disintegrase" to make new attP sites in the Drosophila genome at locations showing chromosomal position effects.

An engineered phiC31 "Disintegrase" able to make an attP site in Drosophila out of an attR-attL pair is described. This was used to generate attP sites at genomic locations where a mini-white (mini-w) transgene was subject to chromosomal position effects (CPE). The first step was random genomic integration of a P-element-based transposon with an insulated mini-w transgene. We then removed the upstream insulator using FLP recombinase to detect CPE. Next mini-w and the downstream insulator were "dis-integrated" leaving behind an attP site. The location is marked by a yellow+ transgene that is flanked by loxP sites, so it can also be removed. Using this system, we generated 10 new attP landing platforms. Three of these showing strong activating CPE were selected for further analysis. We show that the attP sites are functional by integrating in plasmids with attB sites. The CPE is recapitulated and can be blocked by insulators. We show that a dimerized 215 bp fragment of the 500 bp BEAF-dependent scs' insulator containing a high affinity BEAF binding site blocks the CPE, while a monomer of the sequence is less effective. This indicates that two BEAF binding sites make a stronger insulator than a single site. This system could be useful for generating attP sites at prescreened sites for other purposes, such as studying CPE in embryos or other tissues or for use with "trapped" enhancers of interest.

Computational Prediction of Position Effects of Human Chromosome Rearrangements.

Balanced and apparently balanced chromosome abnormalities (BCAs) have long been known to generate disease through position effects, either by altering local networks of gene regulation or positioning genes in architecturally different chromosome domains. Despite these observations, identification of distally affected genes by BCAs is oftentimes neglected, especially when predicted gene disruptions are found elsewhere in the genome. In this unit, we provide detailed instructions on how to run a computational pipeline that identifies relevant candidates of non-coding BCA position effects. This methodology facilitates quick identification of genes potentially involved in disease by non-coding BCAs and other types of rearrangements, and expands on the importance of considering the long-range consequences of genomic lesions.

Regulated Gene Expression as a Tool for Analysis of Heterochromatin Position Effect in Drosophila.

Position effect variegation (PEV) is a perturbation of genes expression resulting from the changes in their chromatin organization due to the abnormal juxtaposition with heterochromatin. The exact molecular mechanisms of PEV remain enigmatic in spite of the long history of PEV studies. Here, we developed a genetic model consisting of PEV-inducing chromosome rearrangement and a reporter gene under control of the UAS regulatory element. Expression of the reporter gene could be regulated by adjustment of the GAL4 transactivator activity. Two UAS-based systems of expression control were tested - with thermosensitive GAL4 repressor GAL80ts and GAL4-based artificial transactivator GeneSwitch. Both systems were able to regulate the expression of the UAS-controlled reporter gene over a wide range, but GAL80ts repressed the reporter gene more efficiently. Measurements of the heterochromatin-mediated repression of the reporter gene in the GAL4+GAL80ts system point to the existence of a threshold level of expression, above which no PEV is observed.

Regulatory variants of FOXG1 in the context of its topological domain organisation.

FOXG1 syndrome is caused by FOXG1 intragenic point mutations, or by long-range position effects (LRPE) of intergenic structural variants. However, the size of the FOXG1 regulatory landscape is uncertain, because the associated topologically associating domain (TAD) in fibroblasts is split into two domains in embryonic stem cells (hESC). Indeed, it has been suggested that the pathogenetic mechanism of deletions that remove the stem-cell-specific TAD boundary may be enhancer adoption due to ectopic activity of enhancer(s) located in the distal hESC-TAD. Herein we map three de novo translocation breakpoints to the proximal regulatory domain of FOXG1. The classical FOXG1 syndrome in these and in other translocation patients, and in a patient with an intergenic deletion that removes the hESC-specific TAD boundary, do not support the hypothesised enhancer adoption as a main contributor to the FOXG1 syndrome. Also, virtual 4 C and HiC-interaction data suggest that the hESC-specific TAD boundary may not be critical for FOXG1 regulation in a majority of human cells and tissues, including brain tissues and a neuronal progenitor cell line. Our data support the importance of a critical regulatory region (SRO) proximal to the hESC-specific TAD boundary. We further narrow this critical region by a deletion distal to the hESC-specific boundary, associated with a milder clinical phenotype. The distance from FOXG1 to the SRO ( > 500 kb) highlight a limitation of ENCODE DNase hypersensitivity data for functional prediction of LRPE. Moreover, the SRO has little overlap with a cluster of frequently associating regions (FIREs) located in the proximal hESC-TAD.

Phenotypic interpretation of complex chromosomal rearrangements informed by nucleotide-level resolution and structural organization of chromatin.

Molecular characterization of balanced chromosomal abnormalities constitutes a powerful tool in understanding the pathogenic mechanisms of complex genetic disorders. Here we report a male with severe global developmental delay in the presence of a complex karyotype and normal microarray and exome studies. The subject, referred to as DGAP294, has two de novo apparently balanced translocations involving chromosomes 1 and 14, and chromosomes 4 and 10, disrupting several different transcripts of adhesion G protein-coupled receptor L2 (ADGRL2) and protocadherin 15 (PCDH15). In addition, a maternally inherited inversion disrupts peptidyl arginine deiminase types 3 and 4 (PADI3 and PADI4) on chromosome 1. None of these gene disruptions explain the patient's phenotype. Using genome regulatory annotations and chromosome conformation data, we predict a position effect ~370 kb upstream of a translocation breakpoint located at 14q12. The position effect involves forkhead box G1 (FOXG1), mutations in which are associated with the congenital form of Rett syndrome and FOXG1 syndrome. We believe the FOXG1 position effect largely accounts for the clinical phenotype in DGAP294, which can be classified as FOXG1 syndrome like. Our findings emphasize the significance of not only analyzing disrupted genes by chromosomal rearrangements, but also evaluating potential long-range position effects in clinical diagnoses.

Variation in Position Effect Variegation Within a Natural Population.

Changes in chromatin state may drive changes in gene expression, and it is of growing interest to understand the population genetic forces that drive differences in chromatin state. Here, we use the phenomenon of position effect variegation (PEV), a well-studied proxy for chromatin state, to survey variation in PEV among a naturally derived population. Further, we explore the genetic architecture of natural variation in factors that modify PEV. While previous mutation screens have identified over 150 suppressors and enhancers of PEV, it remains unknown to what extent allelic variation in these modifiers mediate interindividual variation in PEV. Is natural variation in PEV mediated by segregating genetic variation in known Su(var) and E(var) genes, or is the trait polygenic, with many variants mapping elsewhere in the genome? We designed a dominant mapping study that directly answers this question and suggests that the bulk of the variance in PEV does not map to genes with prior annotated impact to PEV. Instead, we find enrichment of top P-value ranked associations that suggest impact to active promoter and transcription start site proximal regions. This work highlights extensive variation in PEV within a population, and provides a quantitative view of the role naturally segregating autosomal variants play in modifying PEV-a phenomenon that continues to shape our understanding of chromatin state and epigenetics.

Identification of an atypical microdeletion generating the RNF135-SUZ12 chimeric gene and causing a position effect in an NF1 patient with overgrowth.

Neurofibromatosis type I (NF1) microdeletion syndrome, which is present in 4-11% of NF1 patients, is associated with a severe phenotype as it is caused by the deletion of NF1 and other genes in the 17q11.2 region. The variable expressivity of the disease makes it challenging to establish genotype-phenotype correlations, which also affects prognosis and counselling. We here describe a 3-year-old NF1 patient with an atypical deletion and a complex phenotype. The patient showed overgrowth, café au lait spots, inguinal freckling, and neurological abnormalities. The extent of the deletion was determined by means of array comparative genomic hybridisation, and its breakpoints were isolated by means of long-range polymerase chain reaction. Sequence analysis of the deletion junction fragment revealed the occurrence of an Alu-mediated recombination that led to the generation of a chimeric gene consisting of three exons of RNF135 and eleven exons of SUZ12. Interestingly, the deletion shares a common RNF135-centred region with another deletion described in a non-NF1 patient with overgrowth. In comparison with the normal RNF135 allele, the chimeric transcript was 350-fold over-expressed in peripheral blood, and the ADAP2 gene located upstream of RNF135 was also up-regulated. In line with this, the deletion causes the loss of a chromatin TD boundary, which entails the aberrant adoption of distal cis-acting regulatory elements. These findings suggest that RNF135 haploinsufficiency is related to overgrowth in patients with NF1 microdeletion syndrome and, for the first time, strongly indicate a position effect that warrants further genotype-phenotype correlation studies to investigate the possible existence of previously unknown pathogenic mechanisms.

Drosophila Histone Demethylase KDM4A Has Enzymatic and Non-enzymatic Roles in Controlling Heterochromatin Integrity.

Eukaryotic genomes are broadly divided between gene-rich euchromatin and the highly repetitive heterochromatin domain, which is enriched for proteins critical for genome stability and transcriptional silencing. This study shows that Drosophila KDM4A (dKDM4A), previously characterized as a euchromatic histone H3 K36 demethylase and transcriptional regulator, predominantly localizes to heterochromatin and regulates heterochromatin position-effect variegation (PEV), organization of repetitive DNAs, and DNA repair. We demonstrate that dKDM4A demethylase activity is dispensable for PEV. In contrast, dKDM4A enzymatic activity is required to relocate heterochromatic double-strand breaks outside the domain, as well as for organismal survival when DNA repair is compromised. Finally, DNA damage triggers dKDM4A-dependent changes in the levels of H3K56me3, suggesting that dKDM4A demethylates this heterochromatic mark to facilitate repair. We conclude that dKDM4A, in addition to its previously characterized role in euchromatin, utilizes both enzymatic and structural mechanisms to regulate heterochromatin organization and functions.

Computational Prediction of Position Effects of Apparently Balanced Human Chromosomal Rearrangements.

Interpretation of variants of uncertain significance, especially chromosomal rearrangements in non-coding regions of the human genome, remains one of the biggest challenges in modern molecular diagnosis. To improve our understanding and interpretation of such variants, we used high-resolution three-dimensional chromosomal structural data and transcriptional regulatory information to predict position effects and their association with pathogenic phenotypes in 17 subjects with apparently balanced chromosomal abnormalities. We found that the rearrangements predict disruption of long-range chromatin interactions between several enhancers and genes whose annotated clinical features are strongly associated with the subjects' phenotypes. We confirm gene-expression changes for a couple of candidate genes to exemplify the utility of our analysis of position effect. These results highlight the important interplay between chromosomal structure and disease and demonstrate the need to utilize chromatin conformational data for the prediction of position effects in the clinical interpretation of non-coding chromosomal rearrangements.

Replication stress affects the fidelity of nucleosome-mediated epigenetic inheritance.

The fidelity of epigenetic inheritance or, the precision by which epigenetic information is passed along, is an essential parameter for measuring the effectiveness of the process. How the precision of the process is achieved or modulated, however, remains largely elusive. We have performed quantitative measurement of epigenetic fidelity, using position effect variegation (PEV) in Schizosaccharomyces pombe as readout, to explore whether replication perturbation affects nucleosome-mediated epigenetic inheritance. We show that replication stresses, due to either hydroxyurea treatment or various forms of genetic lesions of the replication machinery, reduce the inheritance accuracy of CENP-A/Cnp1 nucleosome positioning within centromere. Mechanistically, we demonstrate that excessive formation of single-stranded DNA, a common molecular abnormality under these conditions, might have correlation with the reduction in fidelity of centromeric chromatin duplication. Furthermore, we show that replication stress broadly changes chromatin structure at various loci in the genome, such as telomere heterochromatin expanding and mating type locus heterochromatin spreading out of the boundaries. Interestingly, the levels of inheritable expanding at sub-telomeric heterochromatin regions are highly variable among independent cell populations. Finally, we show that HU treatment of the multi-cellular organisms C. elegans and D. melanogaster affects epigenetically programmed development and PEV, illustrating the evolutionary conservation of the phenomenon. Replication stress, in addition to its demonstrated role in genetic instability, promotes variable epigenetic instability throughout the epigenome.

Copy number variants analysis in a cohort of isolated and syndromic developmental delay/intellectual disability reveals novel genomic disorders, position effects and candidate disease genes.

BACKGROUND: Array-comparative genomic hybridization (array-CGH) is a widely used technique to detect copy number variants (CNVs) associated with developmental delay/intellectual disability (DD/ID). AIMS: Identification of genomic disorders in DD/ID. MATERIALS AND METHODS: We performed a comprehensive array-CGH investigation of 1,015 consecutive cases with DD/ID and combined literature mining, genetic evidence, evolutionary constraint scores, and functional information in order to assess the pathogenicity of the CNVs. RESULTS: We identified non-benign CNVs in 29% of patients. Amongst the pathogenic variants (11%), detected with a yield consistent with the literature, we found rare genomic disorders and CNVs spanning known disease genes. We further identified and discussed 51 cases with likely pathogenic CNVs spanning novel candidate genes, including genes encoding synaptic components and/or proteins involved in corticogenesis. Additionally, we identified two deletions spanning potential Topological Associated Domain (TAD) boundaries probably affecting the regulatory landscape. DISCUSSION AND CONCLUSION: We show how phenotypic and genetic analyses of array-CGH data allow unraveling complex cases, identifying rare disease genes, and revealing unexpected position effects.

Functional variation in allelic methylomes underscores a strong genetic contribution and reveals novel epigenetic alterations in the human epigenome.

BACKGROUND: The functional impact of genetic variation has been extensively surveyed, revealing that genetic changes correlated to phenotypes lie mostly in non-coding genomic regions. Studies have linked allele-specific genetic changes to gene expression, DNA methylation, and histone marks but these investigations have only been carried out in a limited set of samples. RESULTS: We describe a large-scale coordinated study of allelic and non-allelic effects on DNA methylation, histone mark deposition, and gene expression, detecting the interrelations between epigenetic and functional features at unprecedented resolution. We use information from whole genome and targeted bisulfite sequencing from 910 samples to perform genotype-dependent analyses of allele-specific methylation (ASM) and non-allelic methylation (mQTL). In addition, we introduce a novel genotype-independent test to detect methylation imbalance between chromosomes. Of the ~2.2 million CpGs tested for ASM, mQTL, and genotype-independent effects, we identify ~32% as being genetically regulated (ASM or mQTL) and ~14% as being putatively epigenetically regulated. We also show that epigenetically driven effects are strongly enriched in repressed regions and near transcription start sites, whereas the genetically regulated CpGs are enriched in enhancers. Known imprinted regions are enriched among epigenetically regulated loci, but we also observe several novel genomic regions (e.g., HOX genes) as being epigenetically regulated. Finally, we use our ASM datasets for functional interpretation of disease-associated loci and show the advantage of utilizing naïve T cells for understanding autoimmune diseases. CONCLUSIONS: Our rich catalogue of haploid methylomes across multiple tissues will allow validation of epigenome association studies and exploration of new biological models for allelic exclusion in the human genome.